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1.0 INTRODUCTION ae 


_ Because much of the indirect effects of the M-X project will result from ete 
recreational activities of the construction workers and military personnel, an ei 
analysis of recreational trends is necessary. The addition of a large number of Boel 
people to the area will require USDA Forest Service, National Park, and BLM aoe 
managers to plan for increased use of campgrounds, picnic areas, and wilderness. - 
This planning can be aided by a model capable of identifying areas most likely to be rane 
impacted and predicting the level of use. Such a model is described below. ml 
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2.0 THE INDIRECT EFFECTS MODEL 


The model used in this analysis is a modification of the familiar "gravity 
model" (Cesario, 1969). This model is used to predict the use (in visitor-days per 
year) of all dispersed and developed recreation sites in the region of interest. It 
takes population estimates of the communities in the study area and beyond, and 
uses quantified estimates of their recreational preferences, as well as travel-time 
from each population center to each recreation area, to predict the use of 
recreation areas. 


Each recreation site is inventoried with respect to the availability and quality 
of certain recreational activities: swimming, camping, picnicking, hiking, off-road 
vehicle (ORV) use, boating, water-skiing, and fishing. Snow skiing and hunting were 
not analyzed because these activities do not follow the primary assumption of the 
gravity model: All things being equal, more people will recreate closer to home 
than farther away. Hunting, particularly big game hunting, is controlled largely by 
the awarding of tags often specific to a region. Skiing is a commercial enterprise 
which, at least in Utah, attracts numerous people from out of state and is localized 
in extent. Minor sports, such as white-water canoeing or bird watching, are not 
evaluated because their numbers are lost in the variance. 


Each state was divided into units of approximately 290 sq mi (17 mi X 17 mi) 
for estimation of dispersed recreation. Travel time (in hours) was estimated from 
each population center to each recreation area. Population centers were defined at 
the township or (in Utah) division level and identified by the name of a city or town 
within the population center. Because the average person in Utah spends approxi- 
mately 9.7 days per year in the activities considered here (Utah SCORP, 1980), the 
amount of visitor-days available is easily calculated. A fraction of these people will 
recreate out of state, and numerous visitors travel into the region. These factors 
are included in the model. 


The model was calibrated by recursively comparing predictions to observed use 
figures for more than 500 of the 983 Nevada/Utah recreation sites and adjusting the 
appeal indexes associated with the activity inventories of the sites. The model 
accounted for 99 percent of the variance in use of the sites, with normally 
distributed residuals. For the analysis, the model was run for the 12-year period 
1982 to 1994, using population projections (ETR-37) for baseline with and without 
other projects, with and without M-X for the Proposed Action, and for all eight 
alternatives. 


This model fits observed recreation patterns very well, but calibration for 
areas for which no data exist can only be based on an inventory of, or speculation 
on, available activities. Experience demonstrated that the inventories alone, 
without corroborating use figures, tend to cause overprediction of low use areas and 
underprediction of popular areas. However, the bulk of the recreation use is in 
USFS areas and state and national parks, for which use data are available. The 
model's estimates are good for those areas in general, however, there are areas for 
which predictions are off. 
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2.1 DERIVATION OF THE MODEL 


The model is based upon the gravity model approach (Cesario, 1969; 1973), 
sometimes called the "travel cost" approach (Smith and Kopp, 1980). This approach 
is reviewed by Vickerman (1974). The formulation of the model resembles that of 
Klein (1979). The model can be best understood in its matrix form: 


A = QR + Q (1) 
where 


Vector of use of all recreation sites (visitor-days); 

Matrix of probability of travel from population centers to recreation 
areas; 

Vector of recreation pool from all population centers (visitor-days); 
Vector of out-of-state use of recreation areas (visitor-days) 


pA O> 
ou 


The matrix Q is derived by multiplying each element of a matrix of the 
probability of travel from a population center j to recreation area i by a weighting 
factor: 


meds aie | (2) 


Re} 
H 


exp (-kx, .) (3) 


J 
which is the probability of travel from j to i if all sites i were equally attractive. 
The numerator is a representation of the assumption that use of recreation areas by 
people from a population center j declines exponentially with travel time, x(ij)* 
(Klein, 1979). The constant k(j) is a coefficient determining the slope of the 
travel-time function for various localities (e.g., people in Nevada might be willing 
to travel farther than people in Utah). 


The denominators in equations (2) and (3) normalize the probability matrices so 
that the column sums equal one. This ensures that all the visitors from each 
population center j are distributed among all the areas. 


The term w(i) in equation (2) is the appeal index (weight) of recreation area i. 
Each area is assigned a weight by solution of the matrix equation 
(4) 


W=SN 


where S is a matrix of available activity quality scores for all recreation areas and 
N is a vector of participation rates for the activities considered in this analysis. 


The activities considered are swimming, picnicking, camping, hiking, off-road —— 
vehicle (ORV) use, boating, water skiing, and fishing. Each recreation area is is 

evaluated with respect to these activities. If a given activity is available, then a 
score greater than zero was assigned; otherwise a score of zero was assigned. 
Descriptions of the recreation areas, as well as 1980-use data, were used to assign 


*Xli) =X = - 4 


scores to each area. These scores are multiplied by the fraction of time spent by an 
average citizen of Utah participating in that activity (the vector N). N was 
obtained from the Utah Statewide Comprehensive Outdoor Recreation Plan (SCOR P) 
(1980), which listed the average number of visitor days (1 person/ 12 hours) spent 
enjoying the activities analyzed. These figures were normalized to give the 
vector N (fraction of use row in Table 2.1-1). 


The appeal index expressed algebraically is 


w. = §s.n (5) 


where o is the index of sport (activity). If s(io ) = 0, then no participation in that 
activity is possible on site i. 


The matrix S is, then, an inventory of available resources on a site and an 
index of that site's relative popularity. Further discussion of S is deferred to the 
section on model calibration. 


The vector R in equation (1) is the set of visitor-days available from each 
population center. The average person in Utah spends 9.7 days per year in the 
activities analyzed (Utah SCORP, 1980). The resource for each population center is 
given by Equation (6). 


= a (6) 
9.7m, (1-t;) 


where 1(j) = population of center j, and t is the fraction of that population assumed 
to recreate in areas other than the ones considered. 


Because Nevada SCORP (1977) lists activity preference by activity occasions, 
those data were not usable for this analysis because activity occasions are not 
directly convertible to visitor-days. All populations in Nevada were assumed to 
behave similarly to those of Utah. New Mexico participation rates were used for 
the Texas/New Mexico analysis. 


The algebraic representation of the model is equation (7) 


p, 255 Ng ae 
a. = > SENG YBa 9.71 (1-t 5) + w; (7) on 
Zp, U8; M,) wy 
where the term is w . is out-of-state visitors to site i. Not all areas get significant 
out-of-state visitation, but the so-called Golden Circle of southern Utah (i.e., Zion, 
Bryce Canyon, and Cedar Breaks) does. Out-of-state use was estimated as a 
fraction of in-state ee kg use and added to the solution for a(i): 


Piz Fiotg _ io" 
rasp TEER EFS) 9.7m (lt) } fy (3) 


where f(i) is eae use fraction for area i. Equation (8) took the value of zero 
for most areas. Estimates of f(i) were obtained from USFS, Park Service, and BLM 
district personnel. 
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2.2 IMPLEMENTATION OF THE MODEL 


The states of Nevada and Utah were sectioned into units of approximately 290 
sq mi (17 mi x 17 mi). Each developed recreation area, wilderness, and state and 
national park was located within a grid unit. Dispersed recreation is referenced to 
each grid unit. Travel time from each population center to each grid unit was 
calculated by identifying the best route and dividing the mileage by average speed. 
Average speed was assumed to be 55 mph on paved roads, 45 mph on gravel, and 35 
mph on unimproved roads. Average speed on trails was assumed to be 4 mph (hiking 
speed). The resulting matrix of travel times (X) was very large: 983 sites and 92 
population centers were used in the Nevada/Utah analysis. Thus, the X, P, and Q 
matrices were of rank (983 x 92), or nearly 100,000 elements. 


Because of the necessity to include rural and urban populations, population 
centers were defined on township or (in Utah) division lines. These centers were 
identified by the name of a city or town in the township. 


The population centers included the operating base at Coyote Spring Valley 
and construction camps. There were 18 camps in Nevada and Utah, 15 camps in 
Texas and New Mexico. Other operating bases were not individually identified; heir 
populations were added to the population center to which they were closest. 


U.S. Census (1980) figures were used for the calibration, and HDR population 
projections for trend and high baseline with and without M-X (ETR-37) were used for 
the analysis. The vector A was expanded to a three-dimensional matrix where 
A (alt,y,1)* is the use of recreation areas by alternative, year, and population level. 
Recreation use for the Proposed Action and all eight alternatives was analyzed for 
the years 1982 to 1994, given trend population baseline, trend baseline plus M-X, 
high population baseline, and high baseline plus M-X (ETR-37). For this, A was 
defined as a matrix of rank (983 x 14 x 4). For most analyses, the difference 
between baseline and baseline plus M-X use is the measure of effect. Sometimes 
the percent change was used. 


The model was encoded in Pascal, and run on HDR's Cyber 173 computer, 
requiring nearly 1.1 hours of CPU time to run. Large arrays were defined as file 
types, and stored on disk in core image format. This step allowed the program to fit 
in approximately 45,000 octal words of central memory; over 2 million words of 
central memory would have been required had arrays been stored in core. 


The results of the analysis are presented and discussed in Chapter 4, EIS, 
ETR-18, and in the appendix of this report (Table A-2, Figures A-! through A-26). 


CALIBRATION 


Use-data for 1980 were obtained for more than 500 sites in Nevada and Utah. 
These included USDA Forest Service data on use of all developed recreation sites as 
well as district totals for dispersed areas. Dispersed use figures were allocated to 
appropriate 17 mi square map units on a proportional basis. The Forest Service data 
were obtained from each regional office. 


*Parentheses denote subscript. 
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National Park Service data were available in the form of total visitor-days and 
a breakdown of backcountry use. BLM data were available for a few sites, other 
estimates were obtained from BLM district personnel. The latter types were rough 
estimates. Data on use of non-Forest Service wilderness were especially difficult to 
obtain. These data were estimates from BLM personnel. 


The model was recursively fit to the data. Such parameters as the travel time 
coefficients, k(j), the out of state travel fractions, f(j), and out of state visitors, 
w (j), were adjusted somewhat. For example, Molyneux (personal communication, 
1981) estimated that 60 percent of use of the portions of Dixie National Forest was 
by out-of-state visitors, mostly from Nevada and California. The model allocated 
approximately 30 percent of the users from Nevada, the rest was added by means of 
the out-of-state term (Equation 8). 


The primary method of calibration was to adjust the inventory scores in 
matrix S by comparison of the observed use A(o) to the predicted use A(p). The 
inventory scores of areas for which the predictions were poor were changed until a 
reasonable agreement was reached. The sum of squares due to the model is obtained 
by difference SSM = SST - SSR (Draper and Smith, 1966) where: 


SSR =Z(A ov Ay)” = residual sum of squares 
SST 


r Ay’ (= A,)In = total sum of squares 


The percent of variation explained by the model is: 
R? = 100 x SSM/SST 


The calibration procedure was repeated until R2- 99 percent. The agreement 
of the model with the data is illustrated in Figure 2.2-1. The residuals were 
normally distributed about zero. Comparison of the model to the data is given in 
Table 2.2-1 in Appendix. 


Although the model fit observed data very well, data were available for only 
half of the areas. Therefore, the error of prediction for the rest of the areas is 
unknown. However, since USFS land is the most heavily used, and since data were 
available for all USFS land, the general predictability of the model is probably good. 
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Figure 2.2-1. 
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Comparison ot observed to predicted 
recreation use of Nevada/Utah sites. 


3.0 DISCUSSION 


Hunting and skiing were not analyzed with this model because these activities 
do not obey gravity model assumptions. Hunters tend to go to areas with abundant 
game or are assigned by tag to specific areas. Skiers have to travel to ski areas. 
Commercial recreation was not analyzed because it is presumed to benefit from 
increased use. 


Although the overall fit of the model was excellent, the residuals ranged from 
very small to quite large. Thus, use of a specific area could be significantly 
overpredicted or underpredicted. 


It is possible to extract the number of persons predicted by the model that 
participate in a given activity at a site: 
c Pijrio%s (9) 
ig x( Ys. n °5 % aromel 
5 Pigs io o? 


J 
where 


ais visitor days per year in activity 
o= activity index 
c= fraction of out-of-state visits participating in activity. 


The other symbols are defined above. 


This step is not currently reliable because the use data did not include activity 
breakdowns. Should such data become available, more detailed analysis would be 
possible. 


Certain improvements could be made to increase theoretical accuracy or 
realism, but they were considered to be unjustified at this stage of development. 
The travel-time coefficient, k(j), could vary with activity, giving a set of 
parameters where k (oj) would be the travel coefficient for the sport o and 
population center (j). This step would require the convolution of the several 
probability distributions or, at the minimum, require o probability matrices. This 
could make the model intractable. 


It can be argued that the activity preference vector N used in the model 
should be expanded to account for state-to-state differences and for the expected 
differences between the current population and construction and military personnel. 
Because no such data were available, it was impossible to incorporate those 
differences. However, a study could be performed by letting N take a variety of 
distributions and observing the effect on the output. 


Finally, it can be argued that the assumption of independence of the activities 
is in error. For example, fishing is often associated with camping. While this is 
true, such distinctions would greatly complicate the model or would require a 
network-type model. 


In conclusion, the model predicts the use of developed and dispersed recreation 
resources very well and allows realistic simulation of future use. It does not allow 
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reliable breakdowns of percent activity at a given area. It can be easily calibrated, 
and appears to be highly appropriate for large-scale or regional planning. 
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